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INTRODUCTION 


Several  well-known  polynomials  have  arisen  as  solutions  to  a  class 
of  linear  differential  equations  of  the  second  order.    The  polynomials 
of  Jacobi,  Hermite  and  La guerre  are  classical  examples  that  may  origi- 
nate in  this  manner.    These  polynomials  become  all  the  more  interesting 
due  to  a  number  of  properties  which  they  have  in  common. 

One  of  the  more  interesting  of  these  properties  is  that  of  orthog- 
onality, with  respect  to  a  weight  funotion,  of  a  set  of  polynomial 
solutions  over  a  fundamental  interval.    The  use  of  this  property,  in 
expanding  an  arbitrary  function  in  a  series  of  these  polynomials,  has 
proved  to  be  a  valuable  contribution  to  the  field  of  mathematics. 

In  this  paper  a  class  of  linear  differential  equations  of  the 
fourth  order  will  be  considered.    The  first  step  will  be  the  deter- 
mination of  conditions  whereby  the  differential  equation  yields  a  set 
of  polynomial  solutions.    The  orthogonalization  of  this  solution  set 
and  the  subsequent  formal  expansion  of  an  arbitrary  function  in  a 
series  of  these  orthogonal  polynomial  solutions  are  further  objectives. 
Conditions  for  achieving  these  objectives  will  be  derived.    The  sepa- 
rate cases  for  the  finite,  semi-infinite  and  infinite  intervals  are 
each  discussed  and  the  requirements  for  each  case  are  determined. 

Analogues  of  the  classical  orthogonal  polynomials  of  Jacobi, 
Hermite  and  La guerre  will  be  obtained  in  this  manner. 
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CHAPTER  I 


POLYNOMIAL  SOLUTIONS 


The  Differential  Equation 
Consider  the  differential  equation 

(1.1)       -poO  ^„,v  +  ^(x)^"'*  +  A0t)<fln'  *  Xnttx)^n  =  Oj 

where  the  coefficients  p,  q,  r,  s  and  t,  none  of  which  vanish  identi- 
cally, are  assumed  to  be  real  polynomials  in  the  real  variable  x.  Let 
An  be  a  polynomial  in  n. 

Form  of  Solutions 

Polynomials  of  the  form 

are  desired  as  solutions  of  the  differential  equation  (1.1).  A 
solution  set   yn  ,  n  =  0,1,2,...,  can  thus  be  written  as 

(  %  -  <3L0IX-  ^~an 

(i.3)  <(      =  a01^  <-a(1x  +  ax^ 

•     ■  • 
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The  forcing  of  the  members  of  (1.3)  as  solutions  of  the  differen- 
tial equation  (1.1)  will  impose  certain  restrictions  on  p,  q,  r,  s  and 

t  of  equation  (1.1). 

The  solution  y0  =  a00.    If  y0  =  a00  is  to  be  a  solution  of 

equation  (1.1),  then 

(1.4)  \otLK)40o-o. 

Since  t(x)  i  0  and  a00  i  0,  then  A0  =  0.  This  requirement  imposes  no 
restrictions  on  p,  q,  r,  s  and  t. 

The  solution  yj.  =  a^x  +  an.    If  yi  =  a01x  +  an  is  to  be  a 
solution  of  equation  (1.1),  then 

(1.5)  Atx)40|  +  X,  tu)(aa,x  +  q„)=  o. 

Since  a01  ^  0,  then  Ai  i  0,  as  s(x)  would  vanish  identically  for 
\1  =  0.    Dividing  equation  (1.5)  by  t(x)  yields 

(1.6)  5(x)  --  -  ^Cadl  x  +  aM)  i  o. 

Hence,  S(x)  is  a  polynomial  of  degree  one,  and  equation  (1.1)  becomes 

(1.7)  |3(x)^/v  +<|(x)^r1",  +  aCx)^M"^tCx)S(^x,,4AK,-tCx)^v1  =  0. 

The  solution  y2  =  a02x2  +  aj^x  +  a22«    ^  ?2  "  ^OZ^  +  a12x  +  a22 
is  to  be  a  solution  of  equation  (1.7),  then 

(1.8)  aa(x)a0i+t(^5(x)(iQ61K  +  aa) 

+  Aj.ttxXaoiXVciaX  +  aai)  =  o. 


Dividing  equation  (1.8)  by  t(x)  yields 


(1.9)      R(jcl  - 

■fc  (x) 


2a 


ox 


Hence,  R(x)  is  a  polynomial  of  degree  two,  at  most.    Equation  (1.7) 
now  becomes 

(l.io)    po0<2„,v  +^(x)^y,m4-t(/)R(x)^y,"-Ht(x)5o<)^^'+-Antc^^M=  o. 

The  solution  y3  =  a-oy^  +  *\y£  +  ^zy*-  *  a33»    If  this  solution  is 
to  satisfy  equation  (1.10),  then 

+  t^)5(x)(3a£>3x24  ial3x  +-al3) 

■*  A3  ttx)(awx3+al3X1  +  aa3X'»-G33)  =  °. 

Dividing  equation  (1.11)  by  t(x)  yields 

(1.12)  QCx)  =  3^ 

=  -,-7-  [Rcx)Cuc>3X4aa13)-'-S(x)(3q<J3X14 iol3x  +  a13) 
+  A3  Ca03  x34-a,3X1+alix  +  a33)J  #  o. 

Hence,  Q(x)  is  a  polynomial  of  degree  three,  at  most,  and  equation 
(1.10)  can  now  be  written  as 

(1.13)  p(*)<^,v  +tCx)^cx)^"'  +  t(x)Rlx^nu-t-tCx^  5(x)cjv,'  +A„t(x>jM  =0. 


The  solution  yi±  =  ap/pc4-  +  aj/pc^  +  a£4X2  +  a^yx  +  a/^.    If  this 
solution  is  to  satisfy  equation  (1.13),  then 

(1.14)  l^pOOaoq  +toO<?cx)(i«4a6Hx  H4IS) 
+  t(x)  R.0c)(uaos  xVtc^x  +  iq.  XH) 

4-  tu)sc*XHa0Vx3  +  3al4x*  +  i4avx+a3J 

Dividing  equation  (1.14)  by  t(x)  yields 

(1.15)  POO  =  ^ 

-fc  Cx.) 


+  5Cx)Cc/a0^x3+3a1<fx2+Ici^x  +  a3(A) 

+  Xw  (4oyx  VQ(<f  x3+  aasxl  <-a3^K+asy)  ^  o. 

Hence,  P(x)  is  a  polynomial  of  degree  four,  at  most,  and  equation 
(1.13)  becomes 

Division  by  t(x)  causes  equation  (1.16)  to  assume  the  form 
(1.17)  P(^^n,v  +  ^Cx)^1"+R(x)^"-hSC><^,'  +  An<^0. 


The  solution  y^  =  a^x^  +  ai^x^-l  +•••+  aj^,  k  =  5»6,7,....  In 
order  for  this  solution  to  satisfy  equation  (1.17),  then 

(1.18)  PCx)rK(k-0(k-2Xk-3^KKk-</+--  +  a^a(k-,)K] 

+  Q(x)[K(k-o(k-j)aOKxt<'3^'---^  (.aCk.3)k] 

The  left  member  of  equation  (1.18)  is  thus  a  polynomial  of  degree  k, 
at  most.  By  a  proper  ohoioe  of  the  coefficients  a^,  J  =  0,1,..., k, 
and  A  k,  the  left  member  can  be  made  to  vanish  identically. 

Summary 

The  differential  equation  (1.17)  has  polynomial  solutions  of  the 
form  of  equation  (1.2)  if  P(x)  is  a  polynomial  of  degree  four,  at  most 
Q(x)  is  a  polynomial  of  degree  three,  at  most;  R(x)  is  a  polynomial  of 
degree  two,  at  most;  S(x)  is  a  polynomial  of  degree  one,  and  Xq  =  0. 


CHAPTER  II 

ORTHOGONALITY  OF  THE  SOLUTION  SET 

Derivation  of  Conditions 
In  the  preceding  chapter  conditions  on  the  polynomial  coefficients 
P,  Q,  R  and  S  of  the  differential  equation  (1.17)  have  been  estab- 
lished suoh  that  a  polynomial  solution  set  of  the  form  of  equation 

(1.2)  exists  for  the  differential  equation  (1.17). 

It  is  now  desired  to  derive  a  set  of  conditions  under  which  the 
solution  set  [yn]t  n  =  0,1,2,...,  will  form  an  orthogonal  system,  with 
respect  to  a  weight  function  p(x),  over  the  fundamental  interval  O*,/^). 

Formation  of  the  Basio  Equation 

Consider  the  system  of  equations: 

where  n  ^  m. 

Multiplication  of  equations  (2.1)  and  (2.2)  by  Oym  and  -  Z)yn, 
respectively,  will  yield  the  system: 

(2.3)  />p^V  7>^"'">  VRV'^V5>'^  +  W><^--°, 
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Addition  of  equations  (2.3)  and  (2.k)  gives 

(2.5)  ^pcv^-W)  ^flV^^^^^^VI*^ 

Let  the  following  change  of  variables  be  made, 

(2.6)  a= 

so  that: 

(2.7)  u'  -  V'l"*"^" 

(2.8)  tc"  -  ^•"cjkv,  ^'V'-^''^-^'''^ 

(2.9)  u'-  =  >,vv +  ^'V'-^W'V  -a^'V- 

Also,  let 

(2.10)  ^s  W'">"^ 
so  that 

(2.11)  IT'=  V"  •J"'"  W^*'" 

The  substitution  of  equation  (2.11)  into  equation  (2.9)  gives 

(2.12)  %n%^m-^^n  *  *m-W 
Substitution  of  equation  (2.10)  into  (2.8)  gives 

(2.13)  >",(j»v,  -^^'"^h  r  -^r- 


Equations  (2.13),  (2.12),  (2.7)  and  (2.6)  substituted  into 
equation  (2.5)  give 

(2.1*0  />f(u"l-a»r,)+/aQ^a|,-cr)^^u'  <-flSu  f  ^n'^)^*Smz°j  **mj 
or  more  simply, 

(2.15)  ^fV" +.0G>u."  f/>Ru  Vs u -  V^'V3^*^*"*""^"!*"1"0'  ■"^ 
Equation  (2.15)  is  to  be  known  as  the  basic  equation  for  orthogonality 

Derivative  Form  of  the  Basic  Equation 

Consider  the  following  relationships: 

(2.16)  (^Pu.")'  f  fopVu"  Vp0L"^ 

(2.17)  Jj^V^'&j'l  (/^'-C/^'Ju'  f  C/>*  -tfirtju'l 

(2.18)  [L/ft-^V/iipyJ^ 

Choose  the  condition 

(2.19)  l(pP)'=pQ. 

The  substitution  of  equations  (2.16),  (2.17),  (2.18)  and  (2.19) 
into  the  basic  equation  (2.15)  will  yield 

(2.20)  ^PL."j'-(/3p)'u"f  (.[/><? -(/^W}'-  ^<?)V 

f       ~(p?y"u  *-(p4)'u.' -(ppy<*'  +psu.-2p?v' 
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Equation  (2.20)  can  be  simplified  to 

(2.21)  (z?Pu.")'+  [cpQ-tp?y]u.']'+{cpf(-(pQy+(ppy,W 
+  C?s-(pRy+(pqy,-(ppr]u.  -  a^p^o' 

Equation  (2.21)  is  the  derivative  form  of  the  basic  equation. 

Integration  of  the  Derivative  Form  of  the  Basic  Equation 

Let  equation  (2.21)  be  integrated,  with  respect  to  x,  from  x  =  c< 
to  x  =  ft  .    That  is,  let 

(2.22)  J  (/3Pu."/dK  +  \{[pQ-(pPy]u.']'dK 

'P  (ft 

{[p\i'(pQ)^(pp)"]u]'d%  +  \[ps-(pR)'-h(pc)y,'(pe)"]udx 

ft  fft 


Performing  the  indicated  integration  in  equation  (2.22)  gives 


ft  S 
(2.23)  f/JPu"J    +  {lpQ-(p?yJu'] 


ft 

h[[pR-(pqy  +  (p?)"]u.}  +  [p5-(pR)'+(pQ)"-(p?y"]u.dx 

4 


-  0}  a  *  m. 
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The  Placing  of  Conditions  on  the  Integrated  Basic  Equation 
It  is  desired  to  obtain  the  result 

which  is  the  usual  definition  of  an  orthogonal  system.  Consequently, 
choose  the  conditions: 

(2.24)  />P-0     at  x^*    and  x 

(2.25)  />Q-(fi?)'  =  O       at  j(m     and  x=^ 

(2.26)  -  (/><?)'  +  (p  ?)"  =  o      at  x=*c   and  x-^ 

(2.27)  ^5  -  (piO'  +(/>$)" -(/>*>>'"  5  O. 

These  conditions,  together  with  the  previously  imposed  condition  of 
equation  (2.19)  which  states 

are  imposed  on  equation  (2.23).    The  result  is 


n  ^  Hi. 

Furthermore,  if  the  parameter  values   An  are  all  distinct, 
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Summary 

By  considering  equation  (2.19)  and  its  derivatives,  it  is  possible 
to  simplify  equations  (2.24),  (2.25),  (2.26)  and  (2.27),  and  obtain: 

(i)  p?  =  o    at  -i  and  fa 

(ii)  (p?)'  -  O       at  *  and  fa 

(iii)  AR  ~(pP)n  =0     at       and  fa 

dv)  ^5  =  (p®'-(f>?y\ 

(V)     ^9    =    2.  C^P) 

These  five  conditions  will  henceforth  be  referred  to  as  orthogonality 
conditions  (i),  (ii),  (iii),  (iv)  and  (v). 

It  has  now  been  established  that  under  these  five  orthogonality 
conditions  the  solution  set  [ynj,  n  =  0,1,2,...,  of  the  differential 
equation  (1.17)  will  form  an  orthogonal  system,  with  respect  to  the 
weight  function  ^(x),  over  the  fundamental  interval  (*<,/9). 

This  may  be  more  explicitly  expressed  by  stating  that  if  f  (x) 
is  assumed  to  be  non-negative  in  the  interval  (-r,^),  then  the  set 
[^D^yn],  n  =  0,1,2,...,  forms  an  orthogonal  system  over  the  funda- 
mental interval  (<*,^). 


CHAPTER  III 


EXPANSION  OF  AN  ARBITRARY  FUNCTION  IN  SERIES 


Assume  that  an  arbitrary  function  f (x)  can  be  expanded  as 

(3.1)         f(x)  =  a0^0  +  a,^,  +  ..•  +  a„yM  +  ",  + 


where  the      are  constants  to  be  determined  and  the  yk  are  solutions 
of  equation  (1.17).    The  determination  of  the  constants  ak, 
k  rz  0,1,2,...,  will  enable  f(x)  to  be  formally  expanded  as  a  series 
in  the  polynomial  solutions  of  differential  equation  (1.17). 
Multiply  both  sides  of  equation  (3.1)  by  />yn  to  obtain 


(3.2)  p-f(*)%„  =  a0p<$0y*        a  »pyS      +    p  + 

Integration  of  equation  (3.2),  with  respect  to  x,  over  the  funda- 
mental interval  of  orthogonality  (*,/?)  will  yield 

(fi  (A  (' 

(3.3)  J  />{u)>dx  =  ^ojPw***  ^J/^.^x +■•• 

-C  °*  ol 

ft 


+  a* 


X  f 


+        I  Pl^d*  + 


13 


14 

Equation  (2.29)  may  now  be  utilized  so  that  every  term  of  the 
right  member  of  equation  (3*3) »  with  the  exception  of  the  integrated 
square,  vanishes.  Thus, 


(3.4) 


Each  an,  n  =  0,1,2,...,  can  now  be  uniquely  determined,  subject  only 
to  the  integrability  of  the  expressions  in  equation  (3»4). 

Hence,  an  arbitrary  function  of  x  can  be  formally  expanded  as  a 
series  in  the  solution  set  [ynJ,  n  =  0,1,2,...,  of  the  differential 
equation  (1.17). 


CHAPTER  IV 


SOLUTIONS  IN  THE  FINITE  INTERVAL 

The  Finite  Interval 
The  fundamental  interval  (°Sf/J)  may  conceivably  extend  in  either 
direction,  or  both  directions,  to  infinity.    These  situations  will  be 
discussed  in  Chapters  V  and  VI.    This  chapter,  however,  will  be  devoted 
to  the  consideration  of  the  finite  interval     %ft) ,  where      <  /3. 

The  Weight  Function  A0(x) 

The  satisfaction  of  the  five  orthogonality  conditions,  developed 
in  Chapter  II  and  stated  in  its  summary,  depends  on  the  choice  of  /3(x). 
Conversely,  the  choice  of  ft)(x)  is  vital  only  with  respect  to  the  part 
it  plays  in  the  satisfaction  of  these  orthogonality  conditions. 

A  close  scrutiny  of  these  orthogonality  conditions  indicates  a 
form  of  P(x)  which  is  sufficient  to  accomplish  this  aim. 

From  orthogonality  condition  (i)  it  is  seen  that  the  possibility 
of  /3(x)  vanishing  at  x  =  a  and  at  x  ■  ft  must  be  included. 

Since  orthogonality  condition  (ii)  states  that  (  JP)'  =  p'?  +  ^p'=  0 
at  oL  and  ft  ,  then  the  possibility  that  ft'(x)  vanishes  at  x  =  -x  and  at 
x  ■  ft  must  not  be  excluded. 

From  orthogonality  condition  (iii),  which  states  that 
PR  -  (ftP)"=  ftR  .  p"?  .  Zp'v'~  p?"  =  0  at  oL  and  ft,  it  is  seen  that 
ft  must  be  such  that  ft"(x)  might  vanish  atx=  ={  and  x  =  /3  . 
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Orthogonality  conditions  (iv)  and  (v)  are  identities  which  state, 
respectively,  that  PS  +pn'-p"?  -  3/)"?'-  3/>'P"-  />P'"  and 

pQ.  =  2^'P  +  Zp?' ,    The  polynomial  nature  of  Q  and  S  indicates  that 
p   must  be  a  factor  of  the  right  members  of  both  of  these  identities. 
This  is  feasible  if  p' ,  p"  and  p'"  can  each  be  expressed  as  the  product 
of  p   and  a  rational  function,  where  the  denominator  of  the  rational 
function  divides  the  numerator  of  the  term  in  which  it  appears. 

Henoe,  consider  a  choice  of 

C^.i)  p  =  (x-^Cx-p)) 

where  g  and  h  are  real  numbers.  Successive  differentiations  of  equation 
(4.1)  yield: 

(k  2)     P'  =  ^(x-*)*"V*-^U(x-^(x-/^-'r(x-^ 

-  ( x-w^r x -a) u f ci'S-')^-/3)24-lg|Jv(^-^)(x-/3)-(-^(^-0(x-^)2-  ] 


-  ry.v)Vv.o^('a-'y^-^x-^)3^3^Ua.-.)(x--0^x-yfl)2 
L  (x-*>3  (x-^)J 

4  39U(^i)(x-0Q*(x-,4)  +  K(M-0fh-x)(x-Q3  7 

Consequently,  the  choice  of  p  =  (x-ot)8(X-/3)h,  with  the  proper 
restrictions  on  g  and  h,  will  serve  to  satisfy  the  orthogonality  con- 
ditions . 
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The  Leading  Coefficient  P(x) 

The  determination  of  g  and  h  in  equation  (4.1)  is  contingent  on 
whether  or  not  P(x)  ■  0  has  ot  and/or  p,  as  roots.    A  number  of  possi- 
bilities exists  with  regard  to  the  roots  of  P(x)  =  0.  Investigation 
of  each  of  these  possibilities  will  determine  the  form  of  the  coeffi- 
cients of  the  differential  equation  (1.17). 

The  assumption  that  P(oQ  ^  0.    If  P(x)  =  0  does  not  have  ex.  as 
a  root,  then  application  of  orthogonality  condition  (i)  yields 

(4.5)  p?  =  (x-'OV/^P  =  o 

at  x  =  oL .    Since  ?(<?c)  f  0,  by  assumption,  and  since   <<<y3  ,  then 

(4.6)  >  o. 
Application  of  orthogonality  oondition  (v)  gives 

(4.7)  (x-4)*(X-/3)t,Q  =  zl<*-«)*(*-/Z)U?]' 
or 

(4.8)  Q  =  HA^)",(x-4)",[^U-^)P+  KU-^P  +  ex-oOCx-^P1]. 

The  polynomial  nature  of  Q  requires  that  (x-*)  divide  the  expression 
in  brackets  in  equation  (4.8),  but  if  this  is  true,  then 

(4.9)  ^0^)P<A)  =  O. 

From  the  inequality  (4.6),  g>0.    By  assumption,  P(<x)  /  0  and 
Equation  (4.9)  is  thus  a  contradiction.    Hence,  the  assumption  that 
P(<x)  /  0  is  not  valid. 


18 


The  assumption  that  P(^)  /  0.    This  assumption  leads  to  a  contra- 
diction in  a  manner  similar  to  that  of  the  preceding  section.    This  can 
be  seen  by  interchanging  the  roles  of       and  fi  ,  and  of  g  and  h  in  the 
proof  of  contradiction  of  the  assumption  in  the  preceding  section. 

Consequently,  P(x)  =  0  must  have  both      and  ft  as  simple  roots, 
at  least. 

The  assumption  that  P(x)  =  0  has  qC  as  a  simple  root.    From  the 
preceding  section  it  is  seen  that  P(x)  =  (x-<x)  (x-yS)U(x)  and,  by  the 
assumption  in  this  section,  U(<x)  /  0. 

Application  of  orthogonality  condition  (i)  thus  yields 

(4.10)  /?P  =  (K-o0*+,(x-£)h+,U*O 

at  x  =  cL.    Since  U(<*l)  ^  0,  by  assumption,  and  since    u.<-fi%  then 

The  application  of  orthogonality  condition  (ii)  will  give 

(4.12)  (p?)'    (<s+0Cx-^-/3)^i(J4(U  +  I)(x-c<)^'(X^)UO 

at  x  =  oL,    If  g>0,  then  (/>?)'=  0  at  x  =  aL.    However,  if  -Kg SO, 
then  (g+l)(oC-y9)U(<^)  =  0.    By  inequality  (4.11),  g>-l.    By  assumption, 
U(*0  /  Oand  ot</3.    The  equation  (g+l)(«-/3)U(«:)  =  0  is  thus  a 
contradiction.  Hence, 

(4.13)  c.  >0. 
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Orthogonality  condition  (iii)  states 
(4.1*0      /dlV  (f#y  -  CX-O^X-^  R-[(X-*)'*,0c-i5)k«"  U  J  " 

+-i(5h>o.*»)(x--0*(x^)uu  +  i(^i)c/-^^rx-j)K+,Lj' 
4-  ix-*)a*,fx-4)h4,,u,,J 

-  u  Cu+  ocx-^c)2  u  -  a(  U+  i)(x-0t)a^-^)u' 

at  x  =  Ifg>l,  then  /'R  -  (/3P)"  =  0  at  x  =  However,  if 

0<gH,  then  -g(g+l)(*-,^2U(<*.)  =  0.    By  inequality  (4.13),  g>0.  By 
assumption,  U(<*)  f  0  and  The  equation  -g(g+l)(o(-^)2U(o<:)  «=  0 

is  thus  a  contradiction .  Hence, 

(4.15)  Q  >  '  • 

It  is  now  seen  that  each  of  orthogonality  conditions  (i),  (ii)  and 
(iii)  are  satisfied  if  equation  (4.15)  holds.    The  application  of 
orthogonality  condition  (iv),  however,  will  result  in  a  contradiction 
of  the  assumption  of  this  section.    This  contradiction  is  established 
in  the  work  that  follows. 
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The  application  of  orthogonality  condition  (iv)  will  give 

(4.16)  (x--0*(X->£s  LCx«-c)*(X^)KR]'.f(x-^1■•(x^)k*•o]H, 

f3^(^+/)fx-v)^-7x-/4;Uf'u,^^(u^^.Xx-o0V^)U"'^ 

*  UCU^i)(>-i)(X-^*,(^*ay4.3^^+l)(X-4),»+'(X^W_,U' 
+  3(U+.)Cw^+l(x^)UU%U^)5+V^)k4,0"'Jj 

or 

(4.17)  J  S  (x-^)"1(x-^)"i[2(x"^x"^^+  ^(x"^^!t"^l^^(?<~o(^(x^^ 

-.(.(^(U+^fx-^^CX'/jyu'-S^+iKx^^x-^)^" 
-  kOi+Odi-iVx—O'U  -3U(^i)(x-j.)3(*-/2)u' 

The  polynomial  nature  of  S  requires  that  (x-«0  divide  the  expression 
in  brackets  in  equation  (4.17),  but  if  this  is  true,  then 

(4.18)  -^(^^-O^-^^U)  -O. 
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From  inequality  (4.15),  g  >1.    By  assumption,  U(*)  j  0  and  °c  <  0  , 
Equation  (4.18)  is  thus  a  contradiction.    Hence,  the  assumption  that 
P(x)  =  0  has  ^  as  a  simple  root  is  not  valid. 

The  assumption  that  P(x)  =  0  has    j  as  a  simple  root.  This 
assumption  is  not  valid  either,  as  can  be  seen  by  interchanging  the 
roles  of  oC  and  A  ,  and  of  g  and  h  in  the  proof  of  contradiction  of 
the  assumption  in  the  preceding  section. 

Verification  that  P  =  A(x-^)  (x-i)  ,  where  A  is  a  constant.  The 
contradictions  of  the  assumptions  made  in  the  four  preceding  sections 
leave  only  one  possibility;  P(x)  ■  0  has  both  -y  and  ^  as  double 
roots.    This  is  seen  to  be  the  oase  since  in  Chapter  I  the  coefficient 
P  was  determined  as  a  polynomial  of  degree  four,  at  most. 

Let  the  five  orthogonality  conditions  be  applied  when  A  is 
chosen  as    0  =  (x-*)g(x-3)h  and  P  =  A(x-^)2(x-/s)2. 

Orthogonality  condition  (i)  requires 

(4.19)  /H*-*y^(K-/Ou  +  i  =  o 

at  x  =  X  and  at  x  =  ft  .    Now  A  ^  0,  since  P(x)  =  0  if  A  =  0.  Also, 
since   « ,3  it  is  seen  that 

(4.20)  ^  >-2. 
and 

(4.21)  k>~2.. 
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Application  of  orthogonality  condition  (ii)  yields 

(4.22)  (pvY  =  ^f^o^^^M^1^^^^"^^'^ 

at  x  =      and  at  x  =  ^.    If  g  >-U  then  (/>?)'-  0  at  x  -  If 
-2<g4-l,  then  (g+2)  =  0.    Since   *c<^  and,  by  inequality 

(4.20),  g>-2,  then  (g+2)k-/3>  ^  0.  Hence, 

(4.23)  >  »  ~ ' . 

Also,  if  h>-i,  then  (^P)'«  0  at  x  =  ^ .    If  -2<h<=-l,  then 
(h+2)(4— 0  =  0.    Since  ^c</3   and,  by  inequality  (4.21),  h>-2,  then 
(h+2)(i-«)  i  0.  Hence, 

(4.24)  U  >  -  l . 

By  applying  orthogonality  condition  (v) ,  it  is  seen  that 

(4.25)  70U  $S  i^(x-^*,fA-/«)k*,[($+^(M  +  a+^x--oji, 

or 

(4.26)  Q=r  2A(k-->Cx-/3)LC^+4C)c-/8>4-(,U+i)Cx-^)J. 
When  orthogonality  condition  (iii)  is  applied,  then 

(4.27)  flk~(p?Y  =rv-^)5rx-^)k  R-  ,a  |/«|*i)rj^(x-4srx73)k4*i 


at  x  =  <  and  at  x  =  ^.  If  g  >  0,  then  pR  -  ( p  P)"  =  0  at  x  =  °< . 
If  -l<g^O,  then 

(4.28)  \\u)  -  A  (3+  0^4- i.)(«c-/3^=  o, 

If  h  >  0,  then  pR  -  (p  ?)"  =  0  at  x  =  /3  .    If  -1  <  h  s  0,  then 

(4.29)  R/yd)  -  rtfuOd^i^-oO1^  o. 

Only  orthogonality  condition  (iv)  now  remains  to  be  satisfied, 
this  condition  is  applied,  then 

(4.30)  (X-«dV/3)K5  -  [(x^^x-^Rj'-^U-^S-^^^iJ"' 
or 

(4.3D      5  =  (x-^-Vx-^)-'  U  (x-*)R.  +cx-^Kx-^)R' 

-  3/^fUti)(W+i.^+i.>(x-^)*(x->i)-^fcrfc+»Vk*-Ofx-^3J. 


The  polynomial  nature  of  S  will  require  that  (x-<*)  (x-^)  divide  the 
expression  in  brackets  in  equation  (4.31).    If  (x-<x)  divides  this 
expression,  then 

(4.32)  ^<f*-ti)  R(«0  -  A  =  oj 
or 

(4.33)  ^  (.*i-/6)lti(<<)  -  A  (<d  +  i)(<2±i.)fr-/a)~i]  =  o. 

So,  since      <  p  (  either  g  =  0  or  R(<*}  -  A(g+l)(g+2)(<-/j)2  =  0.  By- 
consideration  of  equation  (4.28)  it  is  seen  that  if  g  =  0,  then 
R(x)  -  A(g+l)(g+2)(*-  ;02  ■  0.  Therefore, 

(4.34)  R(*)=  A^+.Ha+j.K*^*    9*- 1. 

If  (x-,J)  divides  the  expression  in  brackets  in  equation  (4.31), 

then 

(4.35)  l^(A~^f<C/s)-Ak  (u^KW-l)^-*)3  =  oy 
or 

(^•36)  «>lRf^  -A(w+.)(k*-^(/a-.01J  =  o. 

Therefore,  since   *c  <  J ,  either  h  =  0  or  R(;0  -  A(h+l)(h+2)(^-^)2  =  0. 
By  consideration  of  equation  (4.29)  it  is  seen  that  if  h  =  0,  then 
R(£)  -  A(h+l)(h+2)(4--<)2  =  0.  Therefore, 

(4.37)  =  AtU-nKm*^^)2-  k>-l. 
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Equation  (4.31)  can  now  be  written  as 


(4.38)    S  = 


9  R-A^^^'K^^-x)(x-^)x 


4  kR  ~AUU+))(>+>.)(w.V2 


and,  under  the  oonditions  imposed  in  equations  (4. 3*0  and  (4.37).  S  is 
a  first  degree  polynomial. 


serve  as  a  weight  function  in  the  orthogonalization  of  the  solution  set 
£yrJ»  n  =  0,1,2,...,  of  the  differential  equation  (1.17)  over  the 
fundamental  interval  (*,£)  where  both  °L  and  /9  are  finite. 

However,  this  choice  of  /3(x)  will  require  a  number  of  restrictions 
on  the  coefficients  P,  Q,  R  and  S  of  the  differential  equation  (1.17). 
With  regard  to  P  it  has  already  been  determined  that 

Consequently,  with  this  determination  of  P,  equation  (4.26)  gives 

The  coefficient  R  is  of  the  form 


Summary 


The  choice  of  J(x)  -  (x-^)g(x-/9)h,  where  g>-l  and  h  >-l,  will 


where,  from  equations  (4.3*4-)  and  (4.37)  respectively, 
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Equation  (4.38)  gives 


2? 

Examples  of  the  Finite  Interval 

Example  1 

Let  g  =  h  =  0,  oC  *  -1  and   A  «  1.  Then 

p  =  A^oVx-i)1  s  /A  (xH-  «*+  <). 

where 

Re-'")  =  Sr/»  =  B.-Bj.  -f-83j 
fUO  ~-  SA  *   B,  +  B3^ 
so  that  B2  =  0  and       =  8A  -  B^.  Thus, 

R  =  8,Xz-k  (*A  -  B(). 
S  =  M,  -i^A)  X  . 
The  differential  equation  is 

(kW-X*7+\)  >w  + ACSx^frJO'^'VLQ.xVttA-S,)]^'' 

Solutions  of  the  form 
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will  be  obtained. 
Consequently, 


1"> 


The  highest  power  of  x  that  appears  is  n.    The  coefficient  of 

xn  is 

[X*  +-(lB(-i4-A)r>  +aiw(w~0  +  <?Ar\U-iXn-i-HAyi(n-.XM-iV»i-3)JC(on  -O. 
Since  a^  /  0,  then 

=  -n[/\n3  +  2^^  +  (B,-^A>  +B,-|*a]  . 

The  choices  of  A  and  Bj  are  arbitrary,  with  the  exception  that 
B-l  ^  12A  since  S  =  0  if  Bj_  =  12A.    Hence,  choose  A  ■  1  and       ■  14  and 
thus  simplify  An  to 


}  =  _y^<.v»+-i)* 


The  differential  equation  can  now  be  written  as 
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For  each  choioe  of  n,  where  n  -  0,1,2,...,  a  polynomial  solution  of 
degree  n  will  be  obtained.    The  determinations  of  a  few  of  these 
solutions  follow. 

Let  n  =  0.  Then 

so  that 


Thus,  a00  is  arbitrary  and 


Let  n  =  1 .  Then 


so  that 


o-  a0l  x  -  *a„  -  o. 

Thus,  a01  is  arbitrary  and  an  «  o.  Hence, 

^  i  -  CxQl  X  . 


Let  n  =  2.  Then 
so  that 

o-c^x** -  3ia^x  ~(nc(OI  Htau)  =  o. 

Thus,  a02  is  arbitrary,  a12  =  0  and  a22  =  -(l/3)a02.  Hence, 
Let  n  ■  3.  Then 

*1&I3X  +Giaj 

V  -  °> 

so  that 

+  14,3  -^23)  -'^(4&3X3Ki)3xSGA3K+-Cii3)=:  O 


Thus,  ^3  is  arbitrary,  a13  -  0,  a23  =  -(3/5)ao3  and  a33  =  0.  Hence, 

^3  -"^03  (x5-  ^x). 

Let  n  ■  4.  Then 
so  that 

K«NX*-fcVuGos  x1  ^c^xa-xc^}  *  ^x(sa6Sx3  +  3d^xl 

+"  i^i4X+  <X3<J  -4*oo(ciOHXvVG,sY3.t-GMXl  -K.3lf  X  =  O, 

Thus  a04  is  arbitrary,        ■  0,  a24  =  -{S/7)a^t        =  0  and 
a/j4  =  (3/35)aQji|.  Consequently, 
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Further  solutions  yn,  n  ■  5,6,7,...,  may  be  obtained  in  a  similar 
manner. 

These  polynomials  form  an  orthogonal  system  over  the  fundamental 
interval  (-1,1)  and  are  analogues  of  the  Jaoobi  polynomials.  Further- 
more, if  the  aQn,  n  =  0,1,2,...,  are  chosen  in  such  a  manner  that 
yn(l)  =  1,  then  the  preceding  solution  set  is  seen  to  be  precisely  the 
Legendre  polynomials. 


Example  2 

Let  g  =  1,  h  =  0,    x  =  0  and  ^  »  1.  Then 

/3  -  X. 

£  s  B> ,  X*  +-  BAX  +  6 3 

where 

R  (0  -  a  A  =  B,  +  3t  *-  83> 

so  that  B2  =  -(Bi  +  ^A)  and  B3  ■  6k.  Thus, 

R  -  B,X.X-(G(+4A)>C  4  tA. 

5  =  ( aB.-to/^x-  ua.-'to/O. 

The  differential  equation  is 

A(x''-iX3+x1^„,v+/viox3-ifcxl  ^ux) <jv,"'  *-[B,  xMb+^x 
Solutions  of  the  form 

will  be  obtained. 
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Consequently, 


A* 

The  highest  power  of  x  that  appears  is  n.    The  coefficient  of 

xn  is 

Since  a0n     0,  then 

Xh  =  -A  [.3D,  -(oOA  +  6  -HaAO-0(H-2.)  +  A(^-'X»-i3C^-i)] 

=  -h[An3  ^i-A^2-  HB,-i^V4  xB  -44, A  ]. 

The  choices  of  A  and  Bj  are  arbitrary,  with  the  exception  that 
Bj  {  20A  since  S  =  0  if  Bj  ■  20A.    Hence,  choose  A  =  1  and  Bj  =  23  and 
thus  simplify  An  to 

The  differential  equation  can  now  be  written  as 

(  V4-2X3+x1)  3„vV -i-fi^-r^X^x}^'"  +  (13X*-27X-M.)^,." 
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For  each  choice  of  n,  where  n  =  0,1,2,...,  a  polynomial  solution 
of  degree  n  will  be  obtained.    The  determinations  of  a  few  of  these 
solutions  follow. 

Let  n  ■  0.  Then 


so  that 

Thus,  aoo  is  arbitrary  and 
Let  n  =  1.  Then 


so  that 

Thus,  aQi  is  arbitrary  and  a^j  =  -(2/3)aoi«  Consequently, 
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Let  n  =  2.  Then 


so  that 


Thus,  a02  is  arbitrary,  a^2  11  -(6/5)«02  and  a22  =  (3/lO)aQ2»  Hence, 
Let  n  ■  3»  Then 


so  that 

-t-(SX-(.X  3Q63*3  -v- 10,  i  X  ^Gaj) -11  ^(Go5X3^ G„X*+  Ql3X  +  =Oj 

Thus,  a03  is  arbitrary,  a±j  ■  -(12/7)«03»  a23  =  (6/?)a03  31x1 
a33  =  -(^/35)a03«  Consequently, 

Further  solutions  yn,  n  «=  4,5,6, may  be  obtained  in  a  similar 
manner. 


CHAPTER  V 


SOLUTIONS  IN  THE  SEMI-INFINITE  INTERVAL 

The  Semi- Infinite  Interval 
The  case  of  the  fundamental  Interval  (~<  „  )  which  extends  to  infin- 
ity in  one  direction  will  be  discussed  in  this  chapter.    Two  situations 
arise  here,  namely,  the  interval  (v,co)  and  the  interval  (-oofv).  It 
will  be  seen  that  both  situations  are  fundamentally  the  same,  and  con- 
sequently, the  discussion  will  be  based  on  the  interval  (<* ,  oo ) ,  where 
*  is  finite. 

The  Weight  Function  0(x) 

As  in  Chapter  IV,  the  vanishing  of  0(x),  A'(x)  and/?"(x)  at  the 
finite  value  v  must  be  considered  if  /J(x)  is  to  play  its  part  in  the 
satisfaction  of  the  five  orthogonality  conditions. 

This  case,  however,  differs  from  that  of  the  finite  interval  in 
that  the  interval  extends  to  infinity  in  one  direction.    For  the  sake 
of  convenience,  the  condition  that 

will  henceforth  be  denoted  by  the  statement  "f(x)  =  0  at<so". 

The  fact  that  P  is  a  polynomial  prohibits  the  vanishing  of  P  as 
x^oo.    Consequently,  in  orthogonality  condition  (i),  p(x)  must  vanish 
as  x-)oo. 
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Orthogonality  condition  (ii)  will  require  that  p'(x)  vanish  as 
x-*o°.    Orthogonality  condition  (iii)  forces  p"(x)  to  vanish  as  x->o°. 

Consideration  of  the  identities  in  orthogonality  conditions  (iv) 
and  (v)  indicates  that  P  again  must  be  a  factor  of  the  right  members 
of  both  of  these  identities.    Consequently,  if        ^"and/j"'are  each 
expressible  as  the  product  of  J  and  a  rational  function,  where  the 
denominator  of  the  rational  function  divides  the  numerator  of  the  term 
in  which  it  appears,  then  p  will  be  such  that  the  orthogonality 
conditions  are  satisfied. 

Consider  a  choice  of 

(5.1)  p  =  e.-kX(x-oO» 

where  g  and  h  are  both  constants,  h  >0  and  g  is  a  real  number. 
Successive  differentiations  of  equation  (5*1)  will  give: 

(5.2)  ^'r-ke"kXCx^>»+qe"k*(x--6V,5  e-u*(x-^|"~uCx~a0^'l 

(5.3)  p"  =  UV^w^-^ke^^Va^y^e-^w^-1 


(X-*)'*  J 

3  -  ux 


J 


L  ,'x-*)3  J. 

So,  with  the  choice  h?0,  and  the  proper  restrictions  on  g, 
p  =  e"hx(x-^)S  will  be  sufficient  to  satisfy  the  orthogonality 
conditions . 


40 


The  Leading  Coeffioient  P(x) 

hv 

The  function  e       and  its  suocessive  derivatives  do  not  vanish 
at  any  finite  value  x  *  u  •    A  reappraisal  of  the  contradictions 
established  in  Chapter  IV  indicates  the  invalidity  of  P(x)  =  0  having 
u.  as  a  root  of  multiplicity  less  than  two.    Replacement  of  the  factor 
(x-£)h,  in  the  p  (x)  of  Chapter  IV,  by  e_hx  will  effect  these  same 
contradictions.    This  can  be  seen  by  noting  that  neither  (x-/S)h  and 
its  derivatives  nor  e-nx  and  its  derivatives  will  ever  vanish  at  the 
finite  value  x  =  <<  » 

Consequently,  P(x)  =  (x-^)2V(x),  where  V(x)  is  a  polynomial  of 
degree  two,  at  most. 

Application  of  orthogonality  conditions.    With  the  choice  of 
p  =  e-hx(x-  *)g»  where  h>  0,  and  with  P  =  (x-x)2V,  orthogonality 
condition  (v)  gives 

(5.5)  e^*(wV>Q  =  iU^Nx^^v]' 

j 

or 

(5.6)  Q=  UK-»t^[rk(X-^V4.(^4^V+CX.^V,3, 
Let  V  =  CiX2  +  c2x  +  C3.  Then 

(5.7)  Q  3  aC>-ct)[-k^-^)(clxJ+c1x  +  c,^(^-w.Xc>%Vcxx+c3^ 


The  coefficient  Q  is  a  polynomial  of  degree  three,  at  most,  so  the  coef- 
ficient of  x2*  in  the  right  member  of  equation  (5-7)  must  vanish.  Hence, 

(5.8)  -lUC^o, 
but  sinoe  h  ^  0,  then       =  0  and 

(5.9)  V  -  caX 

Application  of  orthogonality  condition  (iv)  yields 

(5.10)  t~UxOt-fl*S  =[e-uMx^)^R]'-[e-UKrx-^+ivJ," 

-  3  Uf«^+-  / )  (<j  f  ^  e  u  *  ( x-     V  -  C  U   4- 1.)  e-  ^  ( x  -  a )  V 1  v 1 
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or 


(5.11)  Ss(X-4*'  [-UX-*)R  i-^R  +CK~^lV+U3(x-oi)3V-3k\<j+0(/~o()?V 

-3U3(x-*)V  +  guC^,x^(>^)v  uL^^M'v1 

Since  S  is  a  polynomial,  then  (x-«)  divides  the  egression  in  brackets 
in  equation  (5.11).    If  (x-<*)  divides  this  expression,  then 

(5.12)  ^  ft 0*^5 C<3  +  .)C<j ■♦-i.')  V(-0-  <3 
and  equation  (5.11)  becomes 

(5.13)  s  s  ^a^a^v  _  UR  ^  +k3(w),v.3u2(a^HA 


From  equation  (5.9),  V  =  C2x  +  C3,  and  thus  v'  ■  C2.  Let 
R  =  Bix2  +  B2x  +  B3,  and  then  equation  (5.13)  becomes 

(5.1ft)    S  5  juy.J^l-^^.K^uCc^^,)  , 

-  3  U5  (x-j)  ^  -i-  3W$+i)0a+^CCi*+c  3}  +  cu(g  t  0  fx -o0cv 


Sinoe  S  is  a  first  degree  polynomial,  the  coefficients  of       and  x2, 

in  the  right  member  of  equation  (5«l2+)»  must  vanish.    If  the  coefficient 

of       vanishes,  then 

(5.15)  k3c1  =  0j 
and  since  h  f  0,  then  C2  =  0,  and  thus 

(5.16)  V  =  c3. 
Also,  if  the  coefficient  of  x2  vanishes,  then 

(5.17)  -UB,  ^U3C3  = 
so  that 

(5.18)  8,=  UaC3. 
Now 


(5.19)  R  -  U7c3Xa  +■  BXX  *  B3) 
and  equation  (5.1*0  can  be  written  as 

(5.20)  s=  ^^^^^Jhia^  -LlB.x+B,) 
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Apply  orthogonality  condition  (i)  to  obtain 

(5.21)  pP  -  c3e~ky  c*~^^  =  o 

at  <*■  and  °o  .    Since  h  >0,  then        ■  0  at  °°  .    For  Z)P  to  vanish  at 
x  =  <^  i  then 

(5.22)  %>"'3l- 
Orthogonality  condition  (ii)  gives 

(5.23)  f^py  =  C3[-Le-^(y^nM^^-u()(^  +  ,J 

Again,  h>0  implies  (  0P)'»  Oat«,  and  if  g>-l,  then  (/3P)'  =  0 
at     .    If  -2<g^-l,  then  g+2  =  0,  which  is  not  possible.  Hence, 

(5.24)  cj  >  _  i . 

The  application  of  orthogonality  condition  (iii)  yields 

(5.25)  OR-         =  e-ky(v^i\ -c3[uVuxcx-^+i 

-  e-kx(x^'[ft.  -c3  kl(*-«Ol+lc3U*+aHx-«0 

at  ^  and  °°  .    Here  h>0  implies  that        -  ( JP)"  =  0  at  oo,  and  if 
g>0,  then  oR  -  (^P)"  =  0  at  <*-.    If  -l<gS0,  then 

(5.26)  Rio^  -  c3(.9  +  ,V^4-ii  =.  o. 


From  equation  (5.12),  if  g  i  0,  then  R(^-)  -  C3(g+l)(g+2)  =  0.  Hence, 

(5-27)  RC-O  =^3C^  +  0C^h-^  j  V1* 

The  Interval  (~3Q,x) 

The  foregoing  development  used  the  fact  that  h  /  0.    If  p  is 
chosen  as  p  =  e_hx(x-^)S,  where  h<  0  (rather  than  h>0),  and  g  >-l, 
then  the  preceding  development  will  be  reproduced  and  will  yield 
comparable  results  for  the  interval  (-  t*>,    .    When  this  interval  is 
under  consideration,  the  condition  that 

will  be  denoted  by  the  statement  "f  (x)  =  0  at-oo». 

Summary 

The  choioe  of    0  ■  e_hx(x-^)g,  where  h>0  (or  h<0),  and  g 
forces  P  to  be  of  the  form 

From  equation  (5*6), 

The  coefficient  R  is  of  the  form 

where,  from  equation  (5.18), 


and,  from  equation  (5»27)> 
From  equation  (5*20), 


These  choices  of  P,  Q,  R  and  S  will  satisfy  all  five  orthogonality 
conditions  for  the  interval  C^,00)  [or  (-^°,^)J  ,  with  respect  to  the 
weight  function  p  =  e~^(x-*)g. 
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Examples  of  the  Semi-infinite  Interval 

Example  1 

Let  h  =  1,  g  =  0,  c^=0  and  consider  the  interval  ).  Then 


Q  :  iJCUCj-C^X*)  r  C  IX1). 

ft  -  B.x1  -f  3>x  ^63 


where 


so  that 


lc3+-  i3x\  +  c3^. 
The  differential  equation  is 

Solutions  of  the  form 


will  be  obtained. 
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Consequently, 

K 

The  highest  power  of  x  that  appears  is  n.    The  coefficient  of 

x11  is 

Since  aQn  /  0,  then 

The  choices  of  C3  and  B2  are  arbitrary,  with  the  exception  that 
B2  /  -^3  since  S  ^  0  if  B2  =  -^3.    Kence,  choose  C3  =  1  and  B2  =  -5 
and  thus  simplify  An  to 

The  differential  equation  can  now  be  written  as 
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For  each  choice  of  n,  where  n  =  0,1,2  a  polynomial  solution 

of  degree  n  will  be  obtained.    The  determinations  of  a  few  of  these 
solutions  follow. 

Let  n  =  0.  Then 


so  that 


Thus,  aQ0  is  arbitrary  and 
Let  n  =  1.  Then 


so  that 


Thus,  a01  is  arbitrary  and  a-^  =  -aQl*  Hence, 


Let  n  =  2.  Then 


UoxYl  K,^  1-  GMj 


so  that 


Thus,  a02  is  arbitrary,  a12  =  -4a02  and  a22  =  2a02.  Consequently, 
Let    Let  n  =  3.  Then 


so  that 


^C-^a3a01xSiQ,3x^a,,)-<i(c<oixVol3>VG1Jx^c<iJ)  =  o 
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Thus,  a03  is  arbitrary,  a^  =  -9ao3>        ~  ^^03  and  a33  =  -^03* 
Consequently, 

u ,  =  a63  fK3- sx2  +ifx-0. 

Further  solutions  yn,  n  =  ^,5»6,...,  can  be  obtained  in  a  similar 
manner. 

The  preceding  set  of  polynomial  solutions  are  orthogonal  over 
the  fundamental  interval  (C,°o),  with  respect  to  the  weight  function 
P  -  e~x,  and  thus  are  considered  as  analogues  to  a  set  of  Laguerre 
polynomials . 
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Example  2 

Let  h  =  -1,  g  =  1,  o<  =1  and  consider  the  interval  (-oe,°<).  Then 

ft*  (*-<)e*. 

Q  s  ltX-i)[c3lx-i)  +  3C3]sC3  (iyJ-  +  xK-v). 

where 

so  that  Bi  =  C3  and  B3  =  5C3  -  B2.  Then 

The  differential  equation  is 

-^C^C3-|3i)]^CT"  +  ^C(j1^3^X  +  C^-MC3)]^(,,  +X^^M  -O. 
Solutions  of  the  form 

will  be  obtained. 


53 


Consequently, 


The  highest  power  of  x  that  appears  is  n.    The  coefficient  of 

xn  is 

Since  aon  f  0,  then 

The  choices  of  C3  and  B2  are  arbitrary,  with  the  exception  that 
B2  f         since  S  h  0  if  B2  =  4C3.    Hence,  choose  C3  =  1  and  B2  =  5 
and  thus  simplify  An  to 

The  differential  equation  can  now  be  written  as 

(y?-:LXM^/VuxV:ix-^)^"'vUl^-x^/'  +  waoj„~o. 
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For  each  choice  of  n,  where  n  =  0,1,2,...,  a  polynomial  solution 
of  degree  n  will  be  obtained.    The  determinations  of  a  few  of  these 
solutions  follow. 

Let  n  =  0.  Then 

so  that 

Thus,  aQ0  is  arbitrary  and 
Let  n  =  1.  Then 


so  that 


Thus,  aQi  is  arbitrary  and  an  =  a.Qi.  Consequently, 
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Let  n  =  2.  Then 
so  that 

Thus,  a02  is  arbitrary,  a^2  a  ^02  and  a22  3  a02'  Consequently, 

=  aoi  (X?  +  H-X+-  [>, 

Let  n  =  3.  Then 

^3  =  aOJx3+a  .^VOijX^a  33, 

so  that 
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Thus,  a03  is  arbitrary,  a  13  =  9^03 »  a23  =  15a-03  and  a33  =  "^03 •  Hence, 

^3  a  ^-03       3  +"c»^a     iSX  -/). 

Further  solutions  yn,  n  =  ^,5»6,...,  can  be  obtained  in  a  similar 
manner. 


CHAPTER  VI 


SOLUTIONS  IN  THE  INFINITE  INTERVAL 

The  Infinite  Interval 
The  fundamental  interval  ipi%0)  may  extend  in  both  directions  to 
infinity.    If  this  is  the  case,  then  a  choice  of  /?(x)  will  have  to 
be  made  that  is  somewhat  different  from  the  /?(x)  of  Chapters  IV  and 
V.    As  in  these  previous  chapters,  the  choice  of  p>(x)  will  determine 
the  basic  forms  of  P,  Q,  R  and  S  of  differential  equation  (1.17). 

The  Weight  Function  p(x) 

The  polynomial  nature  of  the  coefficients  P,  Q,  R  and  S  of  the 
differential  equation  (1.17)  prohibits  their  vanishing  at  ±00.  An 
inspection  of  orthogonality  conditions  (i),  (ii)  and  (iii)  will  then 
show  that  p>%  />'and  p"  must  each  vanish  at  ±°°. 

The  identities  in  orthogonality  conditions  (iv)  and  (v)  will 
again  require  that  p't  p"  and  p'"  be  expressed  as  products  of  f>  and 
a  rational  function,  as  in  the  preceding  cases  for  intervals  with 
at  least  one  finite  end  point.    Also,  the  denominators  of  these 
rational  funotions  must  divide  the  numerators  of  the  terms  in  which 
they  appear. 

Consider  a  choice  of 

(6.1)  jO  =  e-Kx*     K  >  o. 
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Successive  differentiations  of  equation  (6.1)  give: 

(6.2)  p'  =  -iM  e-kx^ 

(6.3)  ^'  ^^l&-^-lU-ux' 

As  the  preceding  equations  indicate,  the  choice  of  p  »  e-^*2, 
h>0,  is  sufficient  to  satisfy  the  orthogonality  conditions. 

Application  of  the  orthogonality  conditions.    In  this  case,  the 
initial  use  of  orthogonality  condition  (iv)  is  revealing.  Since 
is  chosen  as  p  =  e-hx2,  condition  (iv)  states 

(6.5)        e-k*'5  =  C£-k*Ztf)'-(£-kxiP)'" 

=  -)Ue-ulR  4-<srKXV-  t-iUe-^P 
-^U2xVKXV  -*fitxe-kxVVe-kxV")' 


Let  P  =  Ajx*4"  +  A2x^  +  A3X2  +  Aj^x  +  A3  and  let  R  =  B^yp-  +  B2x  +  By 
Solving  equation  (6.5)  for  S  gives 

(6.6)  £  2  -}WX.C6,x%  B^B^  ^  (xO,\-m3^ -iiuVa.X^AxX3 

4.  §u3x3  fA.x^  +  Axx3+  A3y'  -hA^x  ^As) 

^VWx(l*A,X2+(.AxX  **A^-  ClHA.X  -«-c.AxV 

As  S  is  a  polynomial  of  degree  one,  the  coefficients  of  x?,  x^, 
x-5,  x\       and  x2  must  vanish  in  the  right  member  of  equation  (6.6). 
The  coefficient  of  x7  is  Sh^.    Since  h>0,  this  implies 

(6.7)  A,-o. 

The  coefficient  of  x6  is  8h3A2,  and  with  h  >0,  then 

(6.8)  Aa=3. 

The  coefficient  of  x^  is  8h3A--},  and  thus 

(6.9)  A,  =  °- 
The  coefficient  of  x^  is  8h^A^.,  and  so 

(6.10)  A  ^  =  O, 

The  coefficient  of  x3  is  -2hB1  +  8h3A^,  so  that 

(6.11)  B,-HU2Ar. 


The  coefficient  of  x2  is  -2hB2,  thus 

(6.12)  &>."0. 

Consequently,  it  is  seen  that  P  =       and  R  =  ^h^Aijx2  +  B3. 
Equation  (6.6)  can  now  be  written  as 

(6.13)  5  e  -au  83\  -  m^A^X. 
Application  of  orthogonality  condition  (i)  gives 

(6.14)  pp    e~uy?Ay  =  o 
at  4  00, 

Orthogonality  condition  (ii)  states 

(6.15)  ppV  *  -aUA^e-^^o 

Consideration  of  orthogonality  condition  (iii)  shows 

(6.16)  (p?)"*  e~uxiR       !L,e-kv'p  +4U*xVu*Vj 

at  i  oc . 

The  only  remaining  orthogonality  condition  is  (v).    Its  application 
yields 

(6.17)  e-u**Q  -  -HUASK  &~u*] 
and  then 

(6.18)  0}  =  -HU  Aa-X. 
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Summary 

The  choice  of  f>  =  e-k*2,  h>0,  forces  P  to  be  of  the  form 
From  equation  (6.18), 
The  coefficient  R  is  of  the  form 
where,  from  equation  (6.11), 

6,  =  f  t^A*-. 

From  equation  (6.13), 

These  choices  of  P,  Q,  R  and  S  will  satisfy  all  five  orthogonality 

conditions  for  the  interval  (-0°,o°),  with  respect  to  the  weight  function 
p  =  e-hx2. 
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Example  of  the  Infinite  Interval 
Let  h  =  1  and  consider  the  interval         o°) .  Then 

P  =  Aj-. 

5  =   -2.(83+ 2/^3-)  X. 

The  differential  equation  is 

Ay^;v-^AJ-X^,",-HCv^xVB3)^,"-i(83+^i-N;X>,  +  >.^^. 
Solutions  of  the  form 

will  be  obtained. 
Consequently, 


A, 

>"*  = 
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The  highest  power  of  x  that  appears  is  n.    The  coefficient  of 

x11  is 

Since  a0n  f  0,  then 

=  -  m  (4  As  h  -  i  B3  -  J  A  5  ). 

The  choices  of       and  B3  are  arbitrary,  with  the  exception  that 
B3  ^  ~2A5  since  S  =  0  if  B3  =  -2Ay    Hence,  choose  A^  =  1  and  B3  =  -4 
and  thus  simplify   An  to 

K~-  -fin2. 

The  differential  equation  can  now  be  written  as 

For  each  choice  of  n,  where  n  =  0,1,2,...,  a  polynomial  solution 
of  degree  n  will  be  obtained.    The  determinations  of  a  few  of  these 
solutions  follow. 

Let  n  =  0.  Then 

/I  II       , ,   111  \V 

so  that 


0'Ct0o  =  c. 


Thus,  aQ0  is  arbitrary  and 


\o  -  ^co- 


Let  n  =  1.  Then 


so  that 


Thus,  aoi  is  arbitrary  and        =  0.  Consequently, 


Let  n  =  2.  Then 


so  that 
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6i  +■  HX(;a01v  +o,1.^-!(.-,(oe)ix"l+aij.x  <*6ii)*<0 


Thus,  a02  is  arbitrary,  a12  =  0  anc^  a22  =  -(1/2)a02*  Consequently, 
Let  n  =  3.  Then 

^3"  =  k4ft3x  +  aa,3j 

so  that 

-t-la13x  +         *Ua63  x3-f  al3x%a 23  x-t-a3j)-o, 
Thus,  a03  is  arbitrary,  a13  =  0,  a23  =  -(3/2)a03  and  a33  =  0.  Hence, 
Let  n  =  Then 

^  =  ^-CV  X^^CL(yX^  +  6(1iyX2  +  Cl34X  +-Clt,4^ 


so  that 

Thus,  ao/4.  is  arbitrary,  a^  ■  0,  &2h  ■  -3ao4,  a^,  =  0  and  a^  =  (3/*+)ao4. 
Consequently, 

Further  solutions  yn,  n  ■  5»6»7>»»«t  can  be  determined  in  a  similar 
manner. 

The  preceding  solution  set  forms  an  orthogonal  system,  with  respect 
to  the  weight  function  e-*2,  over  the  interval  f00,^).    This  set  of 
polynomials  is  thus  considered  as  an  analogue  to  a  set  of  Hermite 
polynomials . 
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